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ABSTRACT 

The NASA Glenn Research Center (GRC), the Department 
of Energy (DOE), and the Stirling Technology Company (STC) 
are developing Stirling convertors for Stirling Radioisotope 
Power Systems (SRPS) to provide electrical power for future 
NASA deep space missions. STC is developing the 55-We 
Technology Demonstration Convertor (TDC) under contract to 
DOE. Of critical importance to the successful development of 
the Stirling convertor for space power applications is the 
development of a lightweight and highly efficient linear 
alternator. 

This paper presents a 3 -dimensional finite element method 
(FEM) approach for evaluating Stirling convertor linear 
alternators. The model extends a magnetostatic analysis 
previously reported at the 35 th Intersociety Energy Conversion 
Engineering Conference (IECEC) to include the effects of the 
load current. STC’s 55-We linear alternator design was 
selected to validate the model. Spatial plots of magnetic field 
strength (H) are presented in the region of the exciting 
permanent magnets. The margin for permanent magnet 
demagnetization is calculated at the expected magnet operating 
temperature for the near earth environment and for various 
average magnet temperatures. These thermal conditions were 
selected to represent a worst-case condition for the planned 
deep space missions. This paper presents plots that identify 
regions of high H where the potential to alter the magnetic 
moment of the magnets exists. 


INTRODUCTION 

NASA GRC, DOE, and STC are developing a Stirling 
convertor for an advanced radioisotope power system to 
provide spacecraft on-board electric power for NASA deep 
space missions. NASA GRC’s role includes an in-house project 
to provide convertor, component, and materials testing and 
evaluation in support of the overall power system development. 
As a part of this work, NASA GRC has established an in-house 
Stirling Research Laboratory for testing, analyzing, and 
evaluating Stirling machines (Ref [1]). Four 55-We convertors 
(TDC 5, 6, 7, and 8) have been built by STC for NASA and are 
currently under test at GRC. A cross-sectional view of the 
55-We TDC is shown in Figure 1. For more information on the 
TDC, see References [2,3,4]. As another part of this work, 
NASA GRC has been developing Finite Element Analysis 
(FEA) and Finite Element Method (FEM) tools for performing 
various linear alternator thermal and electromagnetic analyses, 
and evaluating design configurations. 

A 3-dimensional (3-D) magnetostatic FEM model of STC’s 
55-We Technology Demonstration Convertor (TDC) linear 
alternator was developed by Geng et al (Ref [5]) to evaluate the 
demagnetization fields affecting the alternator magnets. The 
model, previously reported at the 35 th IECEC, was an open- 
circuit model. The effect of load current flow through the 
alternator windings was not included in the preliminary 
analysis. The predicted alternator open-circuit voltage 
compared well with the experimentally measured value, which 
tended to verify the accuracy of the model. The model was then 
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used to perform a preliminary evaluation of the demagnetization 
fields acting on the magnets. The preliminary evaluation 
showed that for the assumed magnet material, UGIMAX 37B, 
and magnet temperature of 75°C, the magnets margin to 
demagnetization was approximately 70 kA/m. The effect of 
current flow through the alternator windings was expected to 
further reduce this margin to demagnetization for the assumed 
magnet material. 



Figure 1 - Cross-section of STC’s 55-We Technology 
Demonstration Convertor (TDC) 


The 3-D magnetostatic FEM model of STC’s 55-We TDC 
linear alternator has been enhanced to include the effect of 
current flow through the alternator coil windings. The current 
waveform is now a boundary condition for the magnetostatic 
model. The linear alternator current and terminal voltage was 
measured in NASA GRC’s Stirling Research Laboratory for 
TDC 7 at full design stroke. The measured current was used as 
a boundary condition in the new linear alternator model. The 
magnetic properties of the magnets used in the TDCs built for 
NASA are now in the model and offer a better resistance to 
demagnetization than the UGIMAX 37B magnets assumed in 
the original model. The enhanced model was used to generate 
spatial plots of the magnetic field strength (H) in the region of 
the permanent magnets. The margin for magnet 
demagnetization was evaluated at various magnet temperatures. 


OBJECTIVE 

The objective of this work was to develop an analytical tool 
that could aid engineers in the design, development, evaluation 
and understanding of advanced linear alternator designs. 


MAGNETOSTATIC/LOAD CURRENT MODEL 


A 3-D magnetostatic model of STC’s 55-We linear 
alternator design was developed using Ansoft’s Maxwell 3-D 
Field Simulator software (Ref [6]). The magnetostatic model 
consisted of a quarter section of the linear alternator as 
described in Ref [5]. The components included the mover 
laminations, Neodymium-Iron-Boron magnets, stator 
laminations, and copper coils. Magnetostatic flux solutions 
were generated for a series of 13 different mover positions over 
a L2 cm mover stroke. For each mover position, the 

instantaneous current flow through the alternator coil windings 
was applied as a boundary condition. The linear alternator 
current was obtained experimentally by GRC for TDC 7. A 
plot of the experimental data is shown in Figure 2. 
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Figure 2 - Experimental Data of Terminal Voltage, Load 
Current, and Mover Position vs. Time for TDC 7 at Full Stroke 


ELECTRO-MAGNETIC ANALYSIS AT ALTERNATOR 
OPERATING TEMPERATURE UNDER LOAD 

A thermal model of STC’s 55-We linear alternator was 
developed using the ANSYS™ finite element software to 
evaluate the alternator magnet temperatures for potential space 
missions. Many assumptions had to be made in preparing the 
thermal model. For instance, the Stirling convertor pressure 
vessel was assumed to have an unobstructed view to a radiation 
sink temperature of -40°C. This sink temperature is a 
conservative estimate for radiation from a body in space in near 
earth environment. These thermal conditions were selected to 
represent a worst-case condition for the planned deep space 
missions. The emissivity and absorptivity of the radiating 
surfaces were set to 0.91. Surface conductance between 
adjacent parts in contact with each other was set to 0.3 1 W/°C- 
cm 2 . The efficiency of the alternator was assumed to be 85%. 
No convection heat transfer was modeled. Assuming a Stirling 
convertor cold-end temperature of 120°C, the thermal model 
predicted that the average magnet temperature would be 
approximately 75°C. At this temperature, the linear alternator 
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magnets have a remanence B r of 1.213 Tesla and a coercive 
force H c of -900 kA/m. 

Figure 3 shows a plot of the magnetic field strength (H) vs. 
distance along a pair of magnets attached to one stator leg for 
the mover positioned at the end of stroke. This is a mover 
position that leaves one of the magnets almost completely 
exposed. The model predicts that this is the most severe 
magnet loading over the entire stroke. The plot shows a 
maximum localized demagnetization field of roughly 725 kA/m 
on the inside surface of the uncovered magnet. At this magnet 
temperature, demagnetization can occur for demagnetization 
fields larger than 884 kA/m. The margin to demagnetization, 
defined as the difference between H at the knee of the normal 
induction curve and the maximum localized demagnetization 
field, is roughly 159 kA/m. 



0 5 10 15 20 25 30 
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Figure 3 - Magnetic Field Strength |H[ vs. Distance along 
Magnets at 75°C; End of Stroke 


Figure 4 shows a plot of H vs. distance along the magnets 
for the mover positioned at mid-stroke. This mover position 
corresponds to the instant in the cycle of the peak-induced 
voltage in the coils due to magnet flux. The plot shows that the 
maximum demagnetization field is about 690 kA/m. The margin 
to demagnetization is 1 94 kA/m. 

Figure 5 shows a plot of the magnetic field strength (H) vs. 
distance along the magnets for the mover positioned at 
maximum current flow. The current lags the induced voltage by 
approximately 40°. Therefore, the mover position shown in the 
figure is offset from its mean position. The plot shows a 
localized maximum demagnetization field of approximately 
660 kA/m. The magnet loading is less severe in this case in 


comparison with the magnet loading when the mover is at the 
end of stroke. The fringing fields in the magnet region near the 
advancing/receding edge of the mover due to the load current 
appears to be less severe than for the case where the magnet 
was almost completely exposed. The margin to demagnetization 
is roughly 224 kA/m. 



Figure 4 - Magnetic Field Strength |H| vs. Distance along 
Magnets at 75°C; Maximum Induced Voltage (Due to Magnet 
Flux) 
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Figure 5 - Magnetic Field Strength |H| vs. Distance along 
Magnets at 75°C; Maximum Load Current 
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Both Figures 4 and 5 show an interesting characteristic of the 
linear alternator due to the load current. For the portions of the 
magnets covered by the mover, the effect of the load current 
tends to decrease the demagnetization field for one magnet 
while it tends to increase the demagnetization field forthe other 
magnet over a half-cycle. The effect reverses itself during the 
second half of the cycle for each magnet. Even though the load 
current can boost the demagnetization field in the covered 
portion of the magnet, the demagnetization field for the 
uncovered portion of the magnet is substantially larger. Figures 
4 and 5 also show a small kink in the demagnetization curve for 
the inside surface of the magnets adjacent to the trailing edge of 
the mover. This kink might represent a minor fringing effect at 
this location. 


MAGNET LOADING AT VARIOUS TEMPERATURES 

The thermal model of STC’s 55-We linear alternator was 
used to evaluate alternator magnet temperatures for various 
Stirling convertor cold-end temperatures. For Stirling convertor 
cold-end temperatures of 100°C and 80°C, the average magnet 
temperatures were calculated to be approximately 62 °C and 
53°C, respectively. The linear portion of the magnet's normal 
induction curve at each temperature was used in the 
magnetostatic model of the linear alternator to evaluate the 
margin to magnet demagnetization as a function of temperature 
at the most severe magnet load condition (i.e. end of stroke). 
For a magnet temperature of 62°C, the maximum localized 
demagnetization field was calculated at 745 kA/m. At this 
magnet temperature, demagnetization can occur for 
demagnetization fields larger than 1018 kA/m. The margin to 
demagnetization is roughly 273 kA/m. For a magnet 
temperature of 53°C, the maximum localized demagnetization 
field was calculated at 760 kA/m. At this magnet temperature, 
demagnetization can occur for demagnetization fields larger 
than 1108 kA/m. The margin to demagnetization is roughly 
348 kA/m. 

Figure 6 shows a plot of three variables as a function of 
magnet temperature. The solid line represents the knee of the 
normal induction curve based on the magnet manufacturer’s 
data. If the magnets are subjected to demagnetization fields 
above this line, some level of magnet demagnetization may 
occur. The long-dashed line in the figure represents the 
maximum localized demagnetization field affecting the magnets 
as predicted by the magnetostatic/load current model. The 
short-dashed line is the predicted margin to magnet 
demagnetization for TDCs 5, 6, 7, and 8. This plot shows that 
the maximum localized demagnetization JieldLfor the alternator 
magnets decreases as temperature increases but at a lower rate 
than the decrease in the magnet’s ability to resist 
demagnetization. As a result, the margin to demagnetization 
decreases as temperature increases. The vertical lines shown on 


the plot indicate the magnet temperatures calculated for the 
various Stirling convertor cold-end temperatures (Tc). 
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Magnet Temperature (°C) 

Figure 6 - Sensitivity of Resistance to Demagnetization, 
Maximum Localized Demagnetization Field Intensity and 
Margin to Demagnetization as a Function of Temperature 

Table 1 summarizes the margin to demagnetization numbers 
presented in this paper and compares these margins on a 
percentage basis calculated in two different ways. The column 
in the Table labeled “Margin Based on Max Localized Demag. 
Field” computes the percentage by taking the difference 
between the demagnetization field limit and the maximum 
localized demagnetization field and dividing by the 
demagnetization field limit. The column labeled “Margin 
Based on Demag. Field Avg. Over Magnet Vol.” computes tlte 
percentage by taking the difference between the 
demagnetization field limit and the demagnetization field 
averaged over the magnet volume and dividing by the 
demagnetization field limit. The demagnetization field limit 
represents the knee of the normal induction curve at the given 
magnet temperature. 

Table 1 illustrates that the margin to demagnetization can vary 
widely depending upon the method chosen to perform the 
calculation. There is a large difference between the maximum 
localized demagnetization field and the volume averaged 
demagnetization field for a given magnet. A magnet that 
appears safe based on its volume averaged demagnetization 
field may, in fact, be in jeopardy. 
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Table 1 - Margin to Magnet Demagnetization 


Avg. Magnet 
Temp. 
(°C) 

Max. Localized 
Demag. Field 
(kA/m) 

Demag. Field Avg. 
Over Magnet Volume 
(kA/m) 

Demag. Field 
Limit 
(kA/m) 

Margin Based On 
Max. Localized 
Demag. Field 
(%) 

Margin Based On 
Demag. Field Avg. 
Over Magnet Vol. 
(%) 

53 

760 

550 

■ 

31 

50 

62 

745 

539 

■ 

27 

47 

75 

725 

527 

884 

18 

40 


CONCLUSIONS 

The results shown by this preliminary 3-D analysis indicate 
that the highest potential for demagnetization appears to be 
along the inside surface of the uncovered magnets. The 
demagnetization field for the uncovered magnet when the 
mover is at the end of stroke is higher than for when the mover 
is at the position of maximum induced voltage or maximum 
load current. The demagnetization field and the magnet’s 
resistance to demagnetization decrease with increasing 
temperature; however, the magnet’s resistance to 
demagnetization decreases at a higher rate. The magnet’s 
margin to demagnetization decreases as magnet temperature 
increases as expected. Based on the results presented in this 
paper, it appears that the 55-We TDCs built for NASA can 
operate safely at magnet temperatures up to 75°C, provided that 
no abnormal load conditions occur. 
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